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Research Accomplishments 


The main objective of the Joint-Research Interchange NCC2-5149 was to develop computer codes 
for accurate simulation of femtosecond pulse propagation in semiconductor lasers and semicon- 
ductor amplifiers [ 1 ]. The code should take into account all relevant processes such as the interband 
and intraband carrier relaxation mechanisms and the many-body effects arising from the Coulomb 
interaction among charge carriers [2]. This objective was fully accomplished. We made use of a 
previously developed algorithm developed at NASA Ames [3]-[5], The new algorithm was tested 
on several problems of practical importance. One such problem was related to the amplification 
of femtosecond optical pulses in semiconductors. These results were presented in several interna- 
tional conferences over a period of three years. 

With the help of a postdoctoral fellow, we also investigated the origin of instabilities that can 
lead to the formation of femtosecond pulses in different kinds of lasers. We analyzed the oc- 
currence of absolute instabilities in lasers that contain a dispersive host material with third-order 
nonlinearities. Starting from the Maxwell-Bloch equations, we derived general multimode equa- 
tions to distinguish between convective and absolute instabilities. We find that both self-phase 
modulation and intensity-dependent absorption can dramatically affect the absolute stability of 
such lasers. In particular, the self-pulsing threshold (the so-called second laser threshold) can oc- 
cur at few times the first laser threshold even in good-cavity lasers for which no self-pulsing occurs 
in the absence of intensity-dependent absorption. These results were presented in an international 
conference and published in the form of two papers. 
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Research Project 


The Interchange No. NCC2-5149 deals with the emerging technology of photonic (or optoelec- 
tronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown 
together on the same chip. To build such devices and subsystems, one needs to model the entire 
chip. PICs are useful for building components for integrated optical transmitters, integrated optical 
receivers, optical data storage systems, optical interconnects, and optical computers. For example, 
the current commercial rate for optical data transmission is 2.5 gigabits per second, whereas the 
use of shorter pulses to improve optical transmission rates would yield an increase of 400 to 1000 
times. The improved optical data transmitters would be used in telecommunications networks and 
computer local-area networks. Also, these components can be applied to activities in space, such as 
satellite to satellite communications, when the data transmissions are made at optical frequencies. 

The research project consisted of developing accurate computer modeling of electromagnetic 
wave propagation in semiconductors. Such modeling is necessary for the successful development 
of PICs. More specifically, these computer codes would enable the modeling of such devices, 
including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which 
there is femtosecond pulse propagation. Presently, there are no computer codes that could provide 
this modeling. Current codes do not solve the full vector, nonlinear. Maxwell’s equations, which 
are required for these short pulses and also current codes do not solve the semiconductor Bloch 
equations, which are required to accurately describe the material’s interaction with femtosecond 
pulses. The research performed under NCC2-5149 solves the combined Maxwell’s and Bloch’s 
equations. 
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An algorithm has been developed that solves the semiconductor Maxwell-Bloch equations [1], without 
making the standard approximations of a slowly- varying envelope (SVEA) and a rotating- wave (RWA). This 
more exact formulation is applied to simulations of the propagation of ultrashort pulses for which the standard 
approximations reach their limits. This development was motivated by the generation of optical pulses as 
short as 8 fs, which has become possible due to recent progress in ultrafast technology. 

Previously, an algorithm was developed for the Maxwell equations [2,3], without making the SVEA 
for calculations of pulse propagation in nonlinear glasses, which exhibit Kerr-like instantaneous nonlinearities. 
It was found that significant differences can occur when the SVEA is not made in Maxwell equations. 
Specifically, light bullets, of 25 fs duration, were found to be stable [4] with the full Maxwell equations, 
whereas previously, calculations with the nonlinear Schroedinger equation had shown them to be unstable 
[5]. More recently [6], using the algorithm developed in references 2 and 3, calculations showed the formation 
of shock waves on the optical carrier wave. Such results are impossible with the SVEA since the carrier wave 
is eliminated from the calculations. 

In this paper, this new algorithm is applied to studies of ultrafast optical pulse propagation in 
nonlinear semiconductor materials, in which many-body effects due to Coulomb interactions are included. In 
addition to the algorithm for Maxwell’s equations, a new algorithm has been developed for the semiconductor 
Bloch equations that does not make the RWA and the two algorithms have been combined into one for the 
coupled semiconductor Maxwell-Bloch equations. In the Bloch equations the relaxation-time approximation 
[1] has been made for the various intraband scattering processes. 

The Maxwell-Bloch equations for pulse propagation in one spatial dimension are the following. 
Assume that the electric field of a pulse that is propagating along the 2 direction is polarized along the x 
axis and ignore the transverse effects, then the Maxwell equations become 


dD x _ dH v 
dt dz 


dHy dE x 
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where P x is the induced polarization. 

In the case of semiconductors, P x is calculated by using the semiconductor Bloch equations [1], (a 
two band model, one conduction band and one valence band). 
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where n\ and n £ are the occupation probabilities for electron and holes of the wave vector k in the conduction 
and valence bands respectively, tu* = (nj+n£-l) is the population inversion, pi t * and p 2 ,t are the dispersive 
and absorptive components of the dipole moment p*, of the wave vector Ar, pjt = pi t jt 4 *P2 t * an< ^ r 
parameters govern various decay processes. The transition energy fcw* is varied over a sufficiently large range 
to accurately describe the interaction of an ultrashort optical pulse with the semiconductor. n|(<) and n^(t) 




are determined by first computing the chemical potentials and p h {t) from n£(<) and n£(*) respectively 
and then using the formula for a Fermi-Dirac distribution to find the quasi-equilibrium Fermi distributions 
n£(<) and n£(<) 

The generalized Rabi frequency Q k = the effective transition energy hA k , which includes 

the band-gap renormalization, and the induced polarization P T are given respectively by 

nt = i(p^.(0 + S^-.,P*). AA k = hu> k -^2V^(n‘ + n b 1 ), P x (t)=^f Pl , k k 2 dk 

q^k q^k 

Initially an algorithm was developed for the simpler optical Maxwell-Block equations for two-level 
atomic systems [7]. A calculation of self-induced transparency was made for a a 10 fs pulse [8]. Figure 1 
shows the electric field of the pulse at several moments during its propagation inside the medium. The top 
curve on the right side in figure 2 shows the corresponding population inversion at some location as the pulse 
goes by. The top curve on the left side in figure 2 shows the resulting population inversion for self-induced 
transparency when the SVEA and the RWA are made. The remaining curves are a comparison of the two 
methods when the atomic transition frequency is detuned away from the optical carrier frequency. The exact 
method is able to capture the off resonance details that the approximate method is incapable of modeling. 

Next the semiconductor Maxwell-Block equations were solved [8] under the simplifications that the 
Coulomb interaction terms were neglected (the free carrier assumption) and that there were no relaxation 
terms in the equations for the evolution of n\ and n£. Figures 3-6 show colliding pulses, including constructive 
and destructive interference. Figure 7 shows gain curves that were obtained under the free carrier assumption. 

Finally figure 8 shows exciton results. Notice the Is and 2s absorption peaks in the case of the 
2 ps dipole decay time. Here the population inversion was specified at minus one, the generalized Rabi 
frequency was used but the transition energy was not renormalized. The presentation shall include additional 
calculations of propagating and colliding pulses in which all the Coulomb and relaxation terms are included. 

We would like to thank Rolf Binder, Optical Sciences Center, University of Arizona, for his many 
helpful comments concerning algorithm development for the semiconductor Bloch equations. 
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Fig. 1. Two-level atom. Self induced transparency. Fig. 2. Two-level atom. Population inversion. 



Fig. 3. Approaching Pulses Fig. 4. Constructive Interference 



Fig. 5. Destructive Interference 


Fig. 6. Separating Pulses 
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We investigate analytically the occurrence of modulation instability in doped fiber lasers and amplifiers using 
a Maxwell-Bloch description for the dopants and without making the usual parabolic-gain approximation. 
We find a new modulation instability occurring near the Rabi frequency, which is not predicted by the conven- 
tional complex Ginzburg- Landau model. We discuss the implications of this new instability for fiber ampli- 
fiers and lasers and analyze the effects of the saturable host absorption on the laser instabilities Atomic 
detuning is shown to significantly enhance the new modulation instability, in both the normal- and the 
anomalous-dispersion regimes. © 1997 Optical Society of America [S0740-3224(97)04610-9] 


1. INTRODUCTION 

The onset of instabilities in various kinds of lasers is gen- 
erally studied by use of a rate-equation model based on 
the Maxwell-Bloch equations suitable for a two-level 
atomic system. 1 Such a model, often referred to as the 
(detuned) Lorenz-Haken model, includes atomic polariza- 
tion dynamics and has been used extensively over the 
past two decades. Its use leads to the concept of the sec- 
ond laser threshold, defined as the pump level at which 
the continuous-wave (cw) operation of the laser becomes 
unstable through a Hopf bifurcation, resulting in a self- 
pulsing output. At higher pump levels the laser can en- 
ter into a chaotic regime through a period doubling or an- 
other route to chaos. 1,2 

The advent of fiber lasers during the late 1980’s forces 
several changes to this standard model of laser instabili- 
ties, mainly because the optical fiber, acting as a host to 
the dopants, introduces group-velocity dispersion (GVD) 
and self-phase modulation (SPM), both of which must be 
incorporated for a proper description of the onset of insta- 
bilities in fiber lasers. 3 In fact, these two phenomena 
lead to an instability, known as the modulation instability 
(MI), even in an undoped and unpumped optical fiber. 4 
It is therefore reasonable to expect that the presence of 
GVD and SPM in the host fiber would change the nature 
of instabilities in fiber lasers in comparison with other 
kinds of lasers (gas and solid-state lasers) that are well 
described by the standard Lorenz-Haken model. Apart 
from this propagation-based instability, several other ex- 
planations for the observed instabilities in rare-earth- 
doped fiber lasers have been reported. It was shown 
theoretically and experimentally that the existence of ion 
clusters in heavily Er-doped fiber lasers leads to single- 
mode cw or self-pulsing behavior, whereas the same 
model is also applicable to dual-wavelength or bipolarized 
lasers. 5,6 Other theories and experiments on Er-doped fi- 
ber lasers have shown self-pulsing, chaos, and antiphase 
dynamics between the different polarization eigenstates 
of the optical field. 7,8 The explanation for the self- 
pulsing behavior of Nd-doped fiber lasers has been re- 


ported to be driven by the dynamics of the two field- 
polarization eigenstates that depend on the birefringence 
of the fiber. 9,10 In this paper, however, we focus on the 
propagation-driven MI phenomenon, and do not consider 
any field-polarization dynamics. 

In recent years, the MI phenomenon has been investi- 
gated in doped (active) fibers used to make lasers and 
amplifiers. 3,4,11,12 When doing so, one has to consider the 
nonlinear interaction of the dopants with the optical field. 
A natural choice is to model the dopants as a two-level 
system with an atomic polarization dephasing time T 2 
and a population relaxation time T x . By far, the most 
popular model employs the parabolic-gain approximation, 
leading to a complex Ginzburg-Landau (CGL) equation 
for the optical field. 3 One study showed that in erbium- 
doped fiber amplifiers, 11 the threshold for MI is consider- 
ably lowered compared with that for undoped fibers. Re- 
cently, Chen et al . 12 included gain dynamics (governed by 
Ti ) as well as a fast saturable absorber in the model and 
discussed the implications of MI for passively mode- 
locked figure-eight lasers. The full atomic polarization 
dynamics (governed by T 2 ) has, however, been neglected 
so far. 

In this paper we investigate the occurrence of MI be- 
yond the Ginzburg-Landau approximation by considering 
the full T 2 dynamics. We introduce the theoretical 
framework in Section 2 and discuss the consequences for 
amplifiers in Section 3. There, we calculate the steady- 
state solutions and derive a dispersion relation for MI. 
The effect of the population relaxation damping time T j 
and the dipole dephasing time T 2 on MI are studied for 
amplifiers. We find that by cooling the fiber amplifier, 
the bandwidth and the strength of the MI can be greatly 
reduced. In Section 4 we deal with fiber lasers and de- 
rive the dispersion relation for MI at resonance. The role 
of saturable absorption in fiber lasers is investigated, and 
we focus on the possibility of MI occurring in the normal- 
dispersion regime. We find indeed such an instability, 
having its origin in the atomic coherence effects related to 
the atomic polarization dynamics. This new instability is 
found to occur at rather low frequencies (—50 MHz) and 
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may explain the self-starting behavior of mode-locked Nd- 
doped lasers. In Section 5 we discuss the effects of de- 
tuning on the occurrence of MI and discuss the differences 
between normal and anomalous operating regimes. 

2. THEORETICAL FRAMEWORK 

Our starting point is a set of Maxwell -Bloch equations 
that describe the propagation of optical fields in a nonlin- 
ear, dispersive medium doped with two-level atoms (or 
ions). We write the electrical field £f(.x, y, z , t) and the 
induced material polarization y, z , t) as 

X(x, y, z, t) = | xF(x, y)A(z, t) 

x exp[i(/?o 2 - )] + c.c., (1) 

z, t) = \iF{x,y)B(z, t) 

x exp[i(/?(>z - u 0 t)] + c.c., (2) 

where x is the polarization unit vector of the light as- 
sumed to be linearly polarized along the x axis, Fix , y) is 
the fiber-mode profile, and /3 0 is the wave number corre- 
sponding to the carrier frequency w 0 . We assume that 
the field-polarization direction is preserved upon propaga- 
tion and that we are dealing with a polarization- 
preserving single-mode fiber. However, most of the re- 
sults are expected to remain qualitatively valid for 
conventional optical fibers. After substituting Eqs. (1) 
and (2) into Maxwell’s equations, and making the slowly - 
varying-envelope and rotating-wave approximations, we 
obtain the following equations for the slowly varying com- 
plex amplitudes A and B (Ref. 4): 


M 

i 1 ip 2 ^ 2 A 


— 

= - B - - aA 7T 


dz 

2 2 2 dt 2 



+ (6 + i ■y)|A| 2 A, 

(3) 

d B 



d t 

= (i<5- 1 )B - iAg, 
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frequency of the laser spectrum. The main assumptions 
in our model are the homogeneously broadened gain me- 
dium and the neglect of spontaneous emission The 
former is not valid for all doped fibers, but for some types 
of glass hosts it is a reasonable assumption. 4 Since we 
are interested in deterministic instabilities, spontaneous 
emission can be neglected without loss of generality. 

There are two distinct origins of the nonlinear effects in 
Eqs. (3)-(5). The fiber nonlinearity y = n 2 io 0 /cA^ ac- 
counts for SPM effects induced by the host, where n 2 is 
the nonlinear refractive index (units m 2 /W), c is the speed 
of light in vacuum, and A efr is the effective fibercore area. 
For completeness we give the relation between n 2 and the 
nonlinear susceptibility y 13) (units of meters squared per 
volt squared) of the fiber 

3 

n 2 = 9 R e [J' (3 W 0 )], (6) 

4 e 0 n*c 

where n is the background refractive index and e 0 is the 
permittivity of the vacuum. The dopant-induced nonlin- 
ear effects are governed by the saturation power P sal , de- 
fined as 

fi 2 cn eoA eff 

^sal = 2sp q 1 ^ 

2m 7\i 2 

where h is Planck’s constant divided by 2tt and ^ is the 
dipole moment of the atomic transition. Note that Eqs. 
(3)-(5) are written in the frame of reference moving with 
group velocity v g = fi\~ l , which means that t - T 
— fiiz t where T is the time in the rest frame. By doing 
this, we eliminate the term fiiid A/ ST) from the left-hand 
side of Eq. (3). 

We now briefly discuss the relation of Eqs. (3M5) with 
the CGL model. 4,12 When the assumption is made that 
the population relaxation time T j is much longer than all 
other lifetimes, we can approximate the actual gain g by 
its steady-state value g s . This allows Eq. (4) to be ex- 
pressed in the Fourier domain as the well-known 
Lorentzian-shaped nonlinear susceptibility: 


T x = g 0 ~ g + Im(A*B)/P. at , (5) 

where g is the gain realized by pumping the dopants, a is 
the optical loss, is the GVD coefficient of the host fiber, 

6 accounts for saturable host absorption, y is the fiber 
nonlinearity, S ~ (<u 0 - a>)T 2 is the scaled detuning be- 
tween the carrier frequency o> 0 and the atomic resonance 
frequency g 0 is the unsaturated gain, and P Mt is the 
saturation power for the dopants modeled as a homoge- 
neously broadened two-level system. We have written 
Eqs. (3)-(5) in such a way that A has units of >/W, B has 
units of \JWL~ l , and g has units of L _1 , where L is the 
length of either the amplifier or the laser cavity. 

It is important to note that Eqs. (3)— (5) are based on a 
traveling-wave description rather than a standing-wave 
approach that is employed in the conventional rate- 
equation analysis. Since we adopt a traveling-wave ap- 
proach, the optical field A (z, t) in Eqs. (3)— (5) in principle 
can represent a very wide spectrum (or many longitudinal 
modes). The detuning parameter S is thus interpreted as 
the mismatch between the gain peak and the dominant 


A( Aa>) 1 - i(Aa>T 2 + S)' 

where Ao> = w - io 0 is the detuning of the spectral com- 
ponent from the carrier frequency. In the CGL model, 
polarization equation (8) is approximated by a Taylor ex- 
pansion near = 0 up to second order, leading to the 
parabolic-gain approximation that is reasonably accurate 
for small values of T 2 . In the time domain this corre- 
sponds to (generally complex) corrections A(l x and A/3 2 of 
the inverse group velocity and the GVD coefficient /? 2 
(Ref. 3): 


1 1 


WxiS)- 2 S2T 2{1 _ ig)2 , 

( 9 ) 


( 10 ) 


From Eq. (9) we see that only at resonance (S = 0) can 
the resulting pulse propagation equation be written in the 
reference frame moving with the new group velocity (/?j 
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+ A/?,)" 1 , because only at resonance is the correction 
A/?, real. After writing the correction to the GVD as 

A/? 2 (0) = i 4 e ig s T 2 2 , the resulting equations of the 
L/LrL, model become 


aA _ 1 . ^ 1 d 2 A 

— - ~(g-a)A + - (b -i/3 2 ) — 
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Jl dt~ g0 ~8- 
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( 11 ) 


(12) 


the fiT y ’ fT (1) “ d (12) are ^ ood approximations of 
the full model only if 7\ is long enough and T 2 is short 

enough. In this paper we explore the shortcomings of the 
CGL model for realistic fiber lasers and amplifiers and 
find interesting behavior outside the realm of the CGL 
modeh We note that the CGL model is only useful for 
amplifiers for which gain saturation can be neglected- 
otherwise, the gam dispersion b would be z dependent’ 
which seems impractical at best. 

The Maxwell-Bloch equations (3)-(5) can be applied to 
both amplifiers and unidirectional <e.g„ ring) lasers. In 
the case of lasers however, one should, in general, solve a 
complicated boundary-value problem to account for the lo- 
calized losses at the cavity mirrors, a task that requires a 
numerical approach. In this paper we adopt the mean- 
intensity approximation by replacing the localized mirror 
osses with a distributed loss incorporated in the total op- 
tical loss o. In the case of amplifiers, such a mean- 
intensity requirement is not valid: the intensity is 

, ° g y f 2 . d ® pendent Because the steady states are so 
fferent for lasers and amplifiers, a modulation stability 

rs S2f.S£J"“ 
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We consider an amplifier (or absorber) of length L with an 
input power P 0 at z = 0. We first find the time- 
independent (steady-state) solution of Eqs. (3)-(5). For- 
mally, it can be written as 

A„(z) = [P Aiz)] l,2 exp{i<f> s (z)], 


B.(2) = 


A,(z)g s (z) 


6 + i 


g>(z) — g 0 


1 + 


Pa(z) 


P„t< 1 + S 2 ) 


(13) 

(14) 

(15) 




Using the imaginaiy part of Eq. (3), we can write the 
phase profile <p,(z) in terms of the power profile P A (z): 

- rf'fo'P A (z') + JV . , i 6 , 

, he real part ofEq. (3). and using E q . (15), one finds 
the following differential equation for the scaled power 
profile /(z) - P A (z)/[P „ t (1 + s 2 )): P 


This differential equation can be solved, resulting in 
following transcendental equation for P A : 


20P Mt (l + S 2 )z 
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where the coefficients C and C. are given by 
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2C. - c 2 - 1 ± [d - c 2) 2 _ 4(Cj _ C2 )] i/2 i (1 < 


C 1 — 


8o 


20P Mt (l + «5 2 ) 2 ' Cz 20P, at (l + s 2 )' 
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In the absence of saturable absorption (6 = 0), the solu 
ti°n of Eq. (17) is implicitly given by 
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which, in the absence of optical loss (o = 0), reduces tc 
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^ Power profile P A U) is found, the gain profib 
fl3)-(16) P ° nZatl0n profile B ' (z) fo »ow from Eqs 

To study the onset of MI, we follow a standarc 
approach by considering the linear stability of the 
steady-state (cw) solution given above. Considering 
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A(z, t) = [(Pq) 172 + «(*, *) + iv{z, t)] 


Pa(z) 


P o 


1/2 


exp[»>,( 2 )]. 


(23) 


g.(z) 


B(z, 0 = ^ + j [(^o) + PU, O + iq(z, 0] 


Pa(z) 


1/2 


exp[i<p,(z)], 


(24) 


£(*, 0 = [g o + *(*, 0] 1 


Pa(*) l’ 1 

Psatd + <* 2 ). 


(25) 


and linearizing Eqs. (3)-(5) in u, v, p, q , and x, we solve 
the resulting five linear equations in Fourier space by in- 
troducing 


y(z, t) = y 0 exp i dzK(z) - iflt 


I 


y = u, v, p, q y x , (26) 

where y 0 is the initial amplitude, O is the frequency, and 
K(z) is the local wave number of the perturbation. The 
resulting dispersion relation for arbitrary detuning S is 
discussed for fiber lasers in Section 5. The dispersion re- 
lation at resonance (S = 0) is given by 


{[2 iK(z) + g s (z) - 40P A (z)][2iK(z) + *,<*)] 

+ /? 2 2 n 2 [n 2 + sgn(/? 2 )n c 2 (2)]}(i - inr 2 ) 
x [(i - inr,)(i - inr 2 ) + id)} 

~ g,(z)[2iK(z) + g,( 2 )](l - iClT 2 ) 
x [i - inr, - /(*)] -g,(z) 
x [2itf(z) + g s (z) - 46P a (z)] 

x [d - inr,)(i - iciT 2 ) + i(z >] 

+ g s 2 (z)[ 1 - inr, - 7(2 )] = 0, (27) 

where Hz) = P A (z)/P Ml , sgn(ft) - ±1, and n c (z) 
= [4yP,*(z)/|/? 2 |] 1/2 is the critical frequency, i.e., the 
maximum frequency for which MI is found to occur in the 
case of anomalous dispersion in a passive fiber. 4 Before 
we examine the implications of Eq. (27) in various re- 
gimes of parameter space, we note that the Rabi fre- 
quency is somewhat hidden: 


Hiu bi (z) 


Pa(z) 


P MiP 1^2 


1/2 


(28) 


The imaginary part of K(z) determines the local gain 
experienced by the perturbation. It is useful to define 
the total integrated gain at frequency fi as 11 


fi(ft) - 


-2 


/: 


dz Im[/f(n, z)], 


(29) 


where the factor 2 converts h ( fl ) to power gain. MI oc- 
curs whenever the wave number K has a negative imagi- 
nary part. In the case of an amplifier, this means the 


perturbation grows faster than the steady-state power, 
whereas for an absorber it means that the perturbation 
dampens less quickly. Dispersion equation (27) reduces 
to the previously reported ones in the appropriate limits. 
In the absence of saturable host absorption ( 6 = 0), the 
dispersion relation from Ref. 11 is obtained in the limit of 
large T x and short T 2 . 

We now consider the occurrence of MI in various re- 
gimes of parameter space for both amplifiers and absorb- 
ers. Because our model has no restrictions with respect 
to the magnitude of the lifetimes T x and T 2 , we can ex- 
plore MI in regimes where the CGL model has no validity. 
For simplicity, we only consider the local perturbation 
gain because the integration in Eq. (29) can be performed 
analytically in a few limiting cases only. 3 We also ignore 
the possibility of saturable host absorption since two- 
photon absorption is relatively weak in silica fibers, and 
other sources of saturable nonlinearity are rarely present 
in amplifiers. When we discuss MI in lasers, we show 
how even relatively small amounts of saturable host ab- 
sorption can affect the MI drastically. 

We start by investigating the effect of the magnitude of 
the dipole lifetime T 2 . For most fiber amplifiers, T 2 is 
estimated to be near 100 fs, corresponding to a wide gain 
spectrum. Because T x is usually in the range 0.1-10 ms, 
the CGL equation is expected to be a good approximation. 
However, by cooling the fiber, the polarization dephasing 
process can be slowed down substantially, making values 
of T 2 ~ 10 ps readily attainable. 13 

In Fig. 1 we show for various values of T 2 in the range 
0.1-10 ps the MI spectrum for a typical fiber amplifier 
with a 30-dB gain, i.e., exp(g 0 L) = 1000. All other pa- 
rameters are given in the caption. Note that the satura- 
tion power P 8at is inversely proportional to T 2 [Eq. (7)]. 
When Eq. (7) is satisfied for each value of T 2 , the Rabi 
frequency remains a constant for all curves: 

= 1.29 x 10" 3 n c . When the dephasing time T 2 is in- 
creased, two trends are observed. 

First, as can be seen in Fig. 1, increasing T 2 leads to a 
shrinkage of the MI bandwidth, whereas the maximum 



Fig. 1. Modulation instability spectrum for an erbium-doped fi- 
ber amplifier at various values of T 2 (indicated in the figure). 
Parameters are g 0 = 6.91 L" 1 , P 0 = 1 mW, 7\ = 0.1 ms, 0 2 
= -20 ps 2 /L, y=3r 1 L“ 1 , and P Mt = 1 mW when T 2 
= 0.1 ps. For undoped fibers, MI occurs up to 17^/2^ 
= 3.9 GHz 
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Fig. 2. MI spectrum for the amplifier of Fig. 1, for even longer 
dephasing times T 2 (indicated in the figure). In the range 21 
< T 2 < 80 ps, MI is totally quenched. When T 2 approaches 
100 ps, the MI spectrum starts to show a narrow, weak peak 
around the Rabi frequency n iUbl /n r = 0.0013. 

gain within the bandwidth decreases rapidly. When T 2 
= 100 fs ( P A /P 8at = 0.01), MI occurs for frequencies up 
to Q c /2ir = 3.9 GHz, while the peak MI strength is found 
near H = f l c /\j2 . The peak strength is very close to the 
analytical value 2yP A , which is found in the CGL limit. 11 
Already when T 2 = 2 ps (P A /P sat = 0.2), the frequency 
band where MI occurs has shrunk -40%, and near T 2 
— 8 ps (P A = 0.8), MI has almost ceased to occur at 
all. Near T 2 = 21 ps (P A /P sal = 2.1), the MI band van- 
ishes completely. Long before that happens, the strength 
of MI is so weak that it is doubtful whether it can be ob- 
served in a single-pass amplifier. 

Second, during this MI spectrum shrinkage, another 
phenomenon is occurring that is directly caused by the 
two-level system since it involves frequencies close to 
^Rabi > as is shown in Fig. 2. Near the Rabi frequency a 
secondary, weak maximum in MI strength begins to form 
for T 2 > 20 ps. This maximum becomes positive near 
T 2 ~ 80 ps (P a /P m t = 8) and grows with T 2 . When T 2 
is increased further, the MI spectrum slowly returns to its 
original width and strength (out of scale in Fig. 2). Near 
T 2 = 11.5 ns iP A /P aat = 1150), the MI spectrum shows 
again positive MI gain around 1I C , while maintaining a 
narrow (but weak) peak close to il R, bi . At the highly im- 
probable value of T 2 - 1 M s ( P A /P Ml = 10 s ) the MI spec- 
trum is very close to the one at T 2 = 100 fs, and we have 
come full circle. 

Thus we find four regimes of T 2 \ in the first regime 
( 100 fs < T 2 < 21 ps), increasing T 2 leads to a total 
quenching of MI. In the second regime (21 ps < T 2 

< 80 ps), no MI occurs, but the gain around the Rabi fre- 
quency is growing. In the third regime (80ps< T 2 

< 11.5 ns), more and more MI occurs around the Rabi 
frequency, while the gain around 12 r is growing toward a 
positive value again. In the fourth regime, approaching 
the long T 2 limit (1 M s < T 2 < *), the MI spectrum re- 
covers fully to its original (small T 2 ) form. The bound- 
aries between these regimes are, of course, strongly de- 
pendent on the power level P A . For higher power levels 
these boundaries rapidly decrease. 


We further note that the MI band near 12 r is insensiti 
to changes in T x , as long as it is accompanied by a chan 
in the saturation power P sat according to Eq. (7). Ho 
ever, if we keep the saturation power constant up 
changing T x (this can be done by adjusting the dipole n 
ment M decreasing T , leads to a stabilization of t. 
lower frequencies and eventually a reduction of MI alt 
gether. 

We emphasize that the narrow MI peak around n Rab 
so weak that it is questionable whether it can be observe 
in an amplifier. In the case of a laser, however, such 
weak gain may build to a substantial instability ov, 
many round trips, as we discuss in the next sectio 
Since it is not common to use an amplifier in the high 
saturated regime, the emergence of the narrow MI bar. 
near the Rabi frequency is not very practical. Note, hov 
ever, that this narrow MI band near the Rabi frequenc 
does not depend on the sign of /? 2 ; in both normal- an 
anomalous-dispersion regimes, this instability emerges a 
relatively high values of 7Y 
Apart from this new (and for realistic systems, ex 
tremely weak) instability, the full Maxwell-Bloch mode 
agrees with the CGL model qualitatively rather well. O 
course, the quantitative differences become larger as tht 
approximations leading to the CGL model (large T } am 
short T 2 ) become more and more inappropriate. In tht 
next section we find that for lasers the situation can b< 
very different. 


LASERS AT RESONANCE 






Equations (3)-(5) also describe the optical field and the 
gain in a laser, when one assumes that all losses can be 
thought of as being distributed along the cavity. Then, 
the steady-state solution is characterized by a 
z -independent power P 0 and gain g s , and can be written 
as 


A,!*) - (P 0 ) 1/2 expft(p,(2)], (30) 

S - i 

B * iZ) = ~ + Ji A sgs. (31) 


gs go 


1 + 


P«td + s 2 ) 


(32) 


Again, from the real and the imaginary part of Eq. (3), the 
following expressions for the laser power P 0 and the 
phase profile <p,(2) are obtained: 


S, - (a - 26P 0 )( 1 + S 2 ), (33) 


IT = yPo + 


1 j 

2 1 1 + ~ 8 2 



(34) 


Since Eq. (33) is quadratic in P 0 [with use of Eq. (32)], in 
principle, two values for the laser power are found. One 
of these is not physical and corresponds to the antilaser, 
which is characterized by a huge gain and almost zero 
power. 


Similar to the amplifier case, we consider small pertur- 
bations u,v,p,q, and x from the cw state, defined as 
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A(z, t) = [(Po)" 2 + »(*, t ) 

+ iv(z, 0 ]exp[»(p,( 2 )], 


(35) 


S - i 

B(z, l) = , ~ rf tCPo) 1 ' 2 + P(2. 0 

+ iq(z, t)]g 8 exp[i<P s (z)]> 


g(z , f) = [g 0 + *)] 


1 + 


P« tCl + <* 2 )J 


(36) 

(37) 


and linearizing Eqs. (3)-(5) in u, v f p, q, andx, we solve 
the resulting five linear equations in the Fourier space by 
introducing 


y(z , t) = y 0 exp [i(Kz - flf)], 


y = u, v, p, q, x, (38) 

where y 0 denotes the initial amplitude of the perturba- 
tion. Note that because because both laser power P 0 and 
g 5 are z independent, the wave number K is also z inde- 
pendent. At resonance (S = 0), the resulting dispersion 
relation reads 


{[2iK + a - 66P 0 ](2iK + a - 20P O ) 

+ /3 2 2 n 2 [n 2 + sgn(/3 2 )fl c 2 ]}(i - «nr 2 ) 

x [(i - inTiHi - inr 2 ) + j 0 ] 

- g„(2iK + a - 2 OP 0 )(1 - inr 2 ) 

x [1 - inr, - / 0 ] - g,(2iK + a - 6 OP 0 ) 

x [d - inr,)(i - inr 2 ) + / 0 ] 

+ g , 2 (i - mr 1 - / 0 ) = o. (39) 

Here, J 0 = Po f Put< sgntft) = ±1, and 
= (4yP 0 /|/8 2 |) 1/2 is the critical frequency, i.e., the maxi- 
mum frequency for which MI is found in the case of 
anomalous dispersion in a passive fiber. 4 This dispersion 
relation is identical to Eq. (27) when one replaces K , P 0 , 
andg, by their z-dependent counterparts, and Eq. (33) is 
used. 

Dispersion relation (39) reduces to the one previously 
reported by Chen et a/., 12 who employ the CGL model, in 
the appropriate limit. 

Before we proceed with examining the implications of 
Eq. (39) in various regimes of parameter space, we note 
that the Rabi frequency is now given by 

/ P \ 1/2 

nR * bi " (p„ t T,T 2 ) ' <40) 

Merely comparing the relative strengths of the critical 
frequency fl c with the Rabi frequency flR* b i does not pro- 
vide much information about the effect of atomic coher- 
ence on MI. The interaction between the fiber nonlinear- 
ity, the GVD, and the two-level system is much more 
involved. 

In contrast with the amplifier case described in the pre- 
vious section, lasers generally operate in the heavily satu- 
rated regime. This means that the instability near the 
Rabi frequency is now more likely to play a significant 
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Q/2ti fTHz) 



-0.7 -0.5 -0.3 *0.1 

net Ml gain [-2 Im(K)] 

Fig. 3. MI analysis for a figure-eight laser. Solid curves indi- 
cate the results of the full model, while dashed curves show those 
of the CGL model. Top figure shows the net MI gain spectra, 
while the bottom figure shows the corresponding trajectory of the 
eigenvalue K on the complex plane. Parameters are a 
= 0.4 L" 1 , g 0 = 6 L" 1 , /? 2 = -0.09 ps 2 L~\ 0 = 0.1 VT 1 L~\ y 
= 0.008 W _l L" 1 , T 2 = 1.27 ps, T x = 10 8 ps, and P Mt = 10 
mW. 


role. Furthermore, many fiber laser systems, e.g., a 
figure-eight laser, contain an effective saturable absorber 
that causes mode locking. 

We first compare the predictions of Eq. (39) with the 
CGL-based expression. 12 To facilitate comparison, we 
used the same parameters as in Ref. 12. In Fig. 3 we 
show the differences for the case of a figure-eight laser. 
Although the trajectories of K in the complex plane as a 
function of frequency f 1 are quite different for a figure- 
eight laser, the resulting net MI gain spectra agree quite 
well, at least in the central region. The frequency range 
over which positive net MI gain occurs is underestimated 
by 10% by the CGL model. Both models show vanishing 
gain at 100 kHz [indicated by the vertical line at H — 0 in 
Fig. 3(a)], whereas the frequency with highest gain is 
near 200 kHz. The Rabi frequency in this case is 
= 55 MHz, and the^critical frequency is n c = 37 GHz. 
They differ by almost three orders of magnitude. Even 
so, MI occurs for frequencies almost twice as large as fl c . 

The results for the dye-laser parameters are shown in 
Fig. 4. Our model predicts that MI occurs in a narrow 
band near 30 GHz, whereas the CGL model predicts no 
instability at all! The Rabi frequency and the criti- 
cal frequency D c are both close to 24 GHz, which explains 
why the interaction between the fiber and the two-level 
system is so highly nonlinear. 

Here we find the first meaningful qualitative difference 
between the full Maxwell-Bloch model and the CGL 
model. Not surprisingly, the laser power in Fig. 4 is —60 
times the saturation power P Mt , which makes the Rabi 
frequency of the same order as the critical frequency fl c . 
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So for the dye laser of Fig. 4 the interaction between the 
two-level system on the one side and the GVD and the 
SPM on the other side cannot be described within a 
parabolic-gain approximation. 

The pump value g 0 at which the cw state loses its sta- 
bility is often classified as the second threshold, 1 as it an- 
nounces the onset of unstable behavior. Similarly, we 
can identify the MI threshold as the gain above which MI 
occurs. At this threshold, MI occurs only at the fre- 
quency corresponding to the peak gain in the MI spec- 
trum, which can be compared with the frequency with 
which perturbations grow at a Hopf bifurcation. In Figs. 
5 and 6 we show the dependence of the MI threshold as a 
function of d for the case of the figure-eight laser (Fig. 5) 
and the dye laser (Fig. 6). The effect of saturable absorp- 
tion is very dramatic in the case of the figure-eight laser 



IWii (THz) 



Fig. 4. MI analysis for a dye laser. Similar as in Fig. 3, except 
for the parameters: o = 0.1 L' 1 , g 0 = 3 L' 1 , p 2 = -0.09 


L~\ 8 = 0.001 W-> L \ y 

= 10 3 ps, and = 1 m W. 


: 0.008 W-'L- 1 , r 2 = 2.45 ps, T, 



FJk J* MI threshold as a function of saturable absorption 8 for 
the fiber laser of Fig. 3. 



Fig. 6. MI threshold as a function of saturable absorption 0 for 
the dye laser of Fig. 4. 



ii/2n 


Fig. 7. Comparison of net MI gain spectra in the absence of 
saturable absorption (0= 0), for the cases of normal and anoma- 
lous dispersion. Other parameters of the fiber laser are the 
same as m Fig. 3. 


parameters: when 6 is larger than 10‘ 6 W -1 L -1 . MI oc- 
curs immediately after the first (lasing) threshold. A 
qualitatively similar dependence is found for the dye la- 
ser, where the MI threshold gain decreases from — 9 to 
1.5 at 6 = 0.1. This feature explains why a relatively 
weak saturable absorber can lead to passive mode lock- 
ing. Although the presence of a saturable absorber is 
evidently very ui eful for the generation of mode-locked 
pulses, it somewhat obscures our investigation of the in- 
teraction between the two-level system and the fiber non- 
linearity and dispersion. This explains why for both the 
figure-eight laser and the dye laser of Figs. 3 and 4, the 
results hardly change if we consider normal dispersion. 

We therefore examine the case 0 = 0, so that we can 
consider the sole interaction between the fiber nonlinear- 
ity and the two-level system occurring in the absence of 
saturable absorption. The interesting question is 
whether atomic coherence can lead to MI in the normal- 
dispersion regime of the fiber. According to the CGL 
model, this is not possible. 

In Fig. 7 we show, again for a fiber laser, but without 
saturable absorption (ring-cavity instead of figure-eight 
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geometry), the net MI gain spectra for 
= ±0.09 ps 2 L" 1 . For anomalous dispersion, we find an 
approximately 4-GHz-wide MI band centered near 2.5 
GHz and a much narrower and much weaker MI band 
centered near 50 MHz. While the MI band near 2.5 GHz 
vanishes in the case of normal dispersion, the narrow low- 
frequency band survives. So, contrary to what the CGL 
model predicts, the presence of dopants can cause MI in- 
stability in the normal -dispersion regime of the fiber. 
Furthermore, in contrast to the amplifier case (Section 2), 
any positive value of -2 Im(iO should be taken seriously, 
since, in a laser, even the smallest growth of a perturba- 
tion may cause a significant change in the output signal 
after many round trips in the cavity. 

Figure 8 shows the peculiar dependence of this new MI 
at normal dispersion when the population relaxation time 
T j is decreased from 10 ^us to 1.375 /xs. Upon decreasing 
T l , the new MI band initially grows stronger, while shift- 
ing to higher frequencies. Decreasing T x further causes 
the band to weaken and finally to vanish abruptly at 
— 1.35 /xs. 

We stress that our results indicate that a fiber ring la- 
ser, operating in the normal-dispersion regime, may show 
unstable behavior at high pump levels, even in absence of 
additional saturable absorbing mechanisms. At reso- 
nance, the strength of the MI in the normal-dispersion re- 
gime is rather weak, which would imply that the instabil- 
ity needs to build up during many round trips in the ring 
laser. In the next section we discuss the effect of detun- 
ing on the strength and the nature of this new instability. 


5. MODULATION INSTABILITY IN 
DETUNED FIBER LASERS 

With various experimental techniques, e.g., through the 
use of gratings, one can force a fiber laser to operate away 
from the gain peak. In our theory this means that we 
have to deal with the effect of detuning 6. Recall that 
there is no CGL version for the detuned case, as the group 
velocity becomes complex in the parabolic-gain approxi- 
mation. For arbitrary detuning, dispersion relation Eq. 
(39) reads 
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Fig. 8. New MI band at normal dispersion as a function of popu- 
lation relaxation time Tj (indicated). Other parameters of the 
fiber laser are the same as those in Fig. 3. 


where 7 0 = Pq/P^X' Equation (39) is recovered by put- 
ting <5=0 and g s = a - 2 $P 0 . Equation (41) can be 
applied as well for amplifiers by treating P 0 , g s , and K as 
z -dependent quantities. 

We now use Eq. (41) to investigate the effect of detun- 
ing on the MI spectra shown in Fig. 7. The introduction 
of S into the problem makes the situation even more com- 
plex. Instead of only two frequencies, i.e., the critical fre- 
quency D c and the Rabi frequency , the problem 
now is governed by the interaction of three frequencies. 
In Figs. 9 and 10 we show the effect of detuning on the 
bandwidth and the strength of MI for normal and anoma- 
lous dispersion, respectively. Clearly, small detunings 
have a large effect on the occurrence of MI, and the sign of 
the detuning also matters. This spectral asymmetry is 
due to the fiber host nonlinearities. When GVD and 
SPM are absent, Eq. (41) is symmetric in detuning S. In 
the anomalous-dispersion regime (Fig. 10), a small value 
of the detuning connects the two MI bands, one owing to 
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Fig. 9. Effect of detuning on the new MI in the normal- 
dispersion regime. Parameters identical to those in Fig. 7. 



Fig. 10. Similar to Fig. 9, except that the laser now operates in 
the anomalous-dispersion regime. 

the passive fiber MI and the other owing to two-level dy- 
namics. Upon increasing the absolute value of the de- 
tuning, the MI bandwidth and the strength increase loga- 
rithmically. At large detuning (6= 0.1), there is no 
distinction between the normal- and anomalous- 
dispersion case. For large detuning the instability is ap- 
parently dominated by the two-level dynamics. 


6. CONCLUSIONS 

We have analyzed the occurrence of modulation instabil- 
ity (MI) in fiber lasers and amplifiers by considering the 
self-phase modulation, group-velocity dispersion, and the 
saturable host absorption. The gain spectrum has been 
fully considered, in contrast to the parabolic-gain approxi- 
mation employed in the complex Ginzburg-Landau 
(CGL) model. We have derived analytical expressions for 
the MI dispersion A"(H) that naturally reduce to previ- 
ously reported research for both lasers and amplifiers. 

For amplifiers, operating not too heavily saturated and 
in absence of saturable absorption, no qualitative differ- 
ences with the CGL description are found, even in re- 
gimes where the basic approximations of that model are 
violated. Quantitatively, however, the differences can be 


G. H. M. van Tartwijk and G. P. AgrawaJ 

quite substantial. We show that by cooling the fiber am- 
plifier and thereby increasing the dipole dephasing time, 
the occurrence of MI can be quenched. For heavily satu- 
rated amplifiers, we find a new instability located in a 
narrow frequency band around the Rabi frequency. The 
CGL model does not predict such an instability. The 
strength of this new instability is very small, and it is 
questionable that its effect can be detected in a single- 
pass amplifier. 

In lasers, a different picture emerges, since any grow- 
ing perturbation may build up over many round trips 
within the laser cavity. Furthermore, the presence of a 
weakly saturable absorbing mechanism is shown to 
greatly enhance the instability. We compare our results 
with those of Chen et aZ., 12 who used the CGL model to 
investigate MI in a dye laser and a figure-eight laser. 12 
Our results for the figure-eight laser agree rather well, 
whereas we find disagreement for the dye laser, which in 
our model is predicted to have an instability of —30 GHz. 
Further indication that the CGL model should be used 
with caution is given when systems without saturable ab- 
sorption are studied: for a fiber ring laser operating in 

the normal-dispersion regime, a narrow MI band of low 
(—50- MHz) frequencies is found, which is not predicted 
by the CGL model. This may explain the self-starting of 
mode-locked Nd-doped fiber lasers. 

The effect of detuning on the strength and the band- 
width of the new instability can be substantial, since non- 
zero detuning effectively introduces a new frequency into 
the problem. Even for a relatively small detuning, the 
strength and the bandwidth of MI increase logarithmi- 
cally, whereas the difference between normal and anoma- 
lous dispersion becomes smaller. The fiber nonlineari- 
ties cause the MI spectrum to become asymmetric with 
respect to detuning. 
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Absolute Instabilities in Lasers with 
Host-Induced Nonlinearities and Dispersion 

Guido H. M. van Tartwijk and Govind P. Agrawal, Fellow, IEEE 


Abstract — We analyze the occurrence of absolute instabilities 
in lasers that contain a dispersive host material with third- 
order noniinearities. Starting from the Maxweli-Bloch equations, 
we derive general multimode equations to distinguish between 
convective and absolute instabilities. We find that both self-phase 
modulation and intensity-dependent absorption can dramatically 
affect the absolute stability of such lasers. In particular, the self- 
pulsing threshold (the so-called second laser threshold) can occur 
at few times the first laser threshold even in good-cavity lasers for 
which no self-pulsing occurs in the absence of intensity -dependent 
absorption. 

Index Terms — Laser stability, nonlinear optics, optical fiber 
lasers, optical Kerr effect, optical pulse generation, optical prop- 
agation in dispersive media. 

I. Introduction 

A LMOST immediately after the advent of the laser, it 
was recognized that laser output can become unstable, 
resulting in irregular power spikes even at a constant pumping 
level | 1 J. Over the last 30 years or so. laser instabilities 
have been studied extensively both from the fundamental and 
applied viewpoints [2], |3J. The fundamental studies have led 
to the flourishing iield of optical chaos. On the applied side, 
the development of techniques for controlling chaos are being 
used to make lasers tailored for specific applications (high 
power, short pulses, clean far field, etc.). 

Since deterministic chaos is studied in a wide variety 7 of 
disciplines, the understanding of laser instabilities can be 
improved by referring to plasma and fluid instabilities that 
have been studied for a long time. A famous example is 
provided by the Lorenz-Haken equations which are named 
after the fluid dynamic ist Lorenz and the laser theorist Haken 
[3], [4]. In fluid dynamics, instabilities are categorized into 
two types: convective and absolute [5]. Convective instabilities 
are characterized by the growth of localized perturbations 
upon propagation inside a nonlinear medium, while absolute 
instabilities exhibit purely temporal dynamics. Absolute laser 
instabilities have been studied for more than 30 years. The 
Lorenz- Haken equations describe the dynamics of a homo- 
geneously broadened gain medium in a unidirectional ring- 
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cavity. Although rarely stated explicitly, the Lorenz-Haken 
equations can only show absolute instabilities. The fundamen- 
tal concepts such as second laser threshold, self-pulsing. Hopf 
bifurcation, and different routes to chaos are all formulated 
within the context of absolute laser instabilities [3]. 

In the last 15 years or so, new laser systems have been 
designed that are not easily modeled by the Lorenz-Haken 
equations. Examples of such lasers are fiber lasers and solid- 
state (c.g.. Ti : sapphire) lasers, which are capable of producing 
ultrashort optical pulses through passive mode locking while 
operating at a constant pump power. What these lasers have in 
common is that the gain is provided by atoms or ions doped 
inside a host material. As a result, the cavity contains not 
only a gain element but also other nonlinear elements, which 
are responsible for nonlinear processes such as self -phase 
modulation (SPM) and intensity- dependent absorption (IDA) 

[6] . Also, group-velocity dispersion (GVD) of the host medium 
plays a nonncgligible role. Because of the dispersive and 
nonlinear effects, evolution of the optical field over a single 
round trip must be considered, contrary to the Lorenz-Haken 
model in which such effects are ignored. This means that the 
convective nature of any instability must be considered while 
discussing instabilities for such lasers. 

A well-known example of a convective instability occurs 
in nonlinear fiber optics [6]. Optical fibers, without any 
gain element and without any longitudinal resonances (no 
cavity), show a convective instability known as the modulation 
instability. When the power of a CW optical beam becomes 
sufficiently large, the combination of SPM and anomalous 
GVD causes the CW beam to break up spontaneously into 
a pulse train (and eventually into optical solitons) whose rep- 
etition rate depends on the fiber parameters. Mathematically, 
a linear stability analysis shows that perturbations of the form 
exp[— i(f!< — Kz)\ grow exponentially as exp(pz) with a 
growth rate g — — Im(/f) that depends on the frequency 
of perturbation. The repetition rate of the resulting pulses 
corresponds to the frequency D for which the growth rate g 
is maximum. 

Adding gain to the system, e.g., by doping the fiber with 
rare-earth ions and pumping it optically, can affect consider- 
ably the conditions under which modulation instability arises 

[7] , The instability, however, remains convective in nature. 
When such a host material (with or without gain) is put into a 
cavity, the resulting boundary conditions at the cavity mirrors 
can change the nature of the instability from convective to 
absolute. Feedback is a necessary ingredient for absolute 
instabilities to occur. A well-known example is the Ikeda 
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instability [8], which arises when a Kerr medium is placed in 
a unidirectional ring cavity. Even without gain and dispersion, 
the feedback mechanism provided by the cavity results in an 
absolute instability. 

In this paper, we discuss under what conditions a convective 
instability becomes absolute in a laser. In Section IL we 
derive, starting from the Maxwell-Bloch equations, a set of 
Lorenz-Ha ken -type multimode equations capable of describ- 
ing the temporal evolution of a laser whose cavity includes 
optical elements exhibiting dispersion and nonlinearities. The 
usefulness of this new set of equations is illustrated in Section 
III by considering a relatively simple case of a single-mode 
laser. We discuss the stability of that mode as a function of 
host nonlinearities. 


II. Maxwell-Bloch Equations 

For definiteness, we focus on a fiber laser although the 
analysis can be applied to any solid-state laser with some 
modifications. Our starting point is a set of Maxwell-Bloch 
equations describing the propagation of optical fields in an 
optical fiber, doped with rare-earth ions. We write the optical 
field £ and the dopant-induced polarization V as 

y, z, t) - la-F(a-, j/). 4(2, <) exp[i(Ao2 - w 0 <)] + c.c. 

( 1 ) 

T(r, y, z , () — y)B(z , t) exp[i(/(o2 - w 0 <)] + c - c - 

( 2 ) 


where r is the polarization unit vector of light assumed 
to be linearly polarized along the x axis. F(x, y) is the 
transverse profile of the fundamental fiber mode, and Kq is 
the wavenumber corresponding to the optical frequency wo- 
We assume that the field-polarization direction is preserved 
upon propagation. After substituting (1) and (2) in Maxwell's 
equations, modeling dopants as a homogeneously broadened 
two-level system, and making use of the slowly varying 
envelope and rota ting-wave approximations, we obtain the fol- 
lowing equations for the slowly varying complex amplitudes 
A and B [6]: 


dA 

dz 


Vy dt 

_ dB 
— 

dt 

T dg 

Tl m 


- t 2 b - \ - l -y ^ + il\A\ 2 A (3) 

= - (1 -i6)B -iAg (4) 

= g D -g + lm{A*B)/P. (5) 


where g is the gain realized by pumping the dopants, a is the 
optical loss of the host fiber, 7\ is the population lifetime 
of the dopants. Tn is the dipole -dephasing time, t ) g is the 
group velocity, /£» is the GVT) coefficient of the host fiber, 
the complex parameter 7 accounts for the host nonlinearities 
responsible for SPM and IDA. 6 = (u>c — 2 is the scaled 

detuning between the optical frequency u^o and the atomic 
resonance frequency u?^, gc is the unsaturated gain, and P t is 
the saturation power for the dopants. We have written (3)-{5) 
in such a way that A has units of \/W* B has units of 
\AV * m” 1 . and g has units of m“ l . 


The main assumptions in our model are the use of a 
homogeneously broadened gain medium and the neglect of 
the stochastic nature of spontaneous emission. The former is 
not valid for all dopants but is a reasonable assumption for 
many types of dopants [6], The latter can be justified if one is 
interested only in deterministic instabilities. 

There are two distinct origins of the nonlinear effects in 
(3)-(5). The host nonlinearity 7 = 7' + *7" accounts for SPM 
and IDA effects induced by the silica fiber. The SPM effects 
arc governed by 7' = n 2 u>o/cAeff, where n 2 is the nonlinear- 
index coefficient, c is the speed of light in vacuum, and A^r is 
the effective mode area [6]. The effects of EDA are accounted 
for by 7". When 7* > 0, the loss in the cavity increases with 
intensity, modeling processes such as two-photon absorption 
[6], In contrast, negative values for 7" imply a decrease in 
cavity losses with increasing intensity and model fast saturable 
absorption. The dopant-induced nonlinear effects are governed 
by the saturation power P t = fi 2 cneo^«ff/(2fi 2 7i72). where 
h is Planck's constant divided by 2x% fi is the dipole moment 
of the atomic transition, and n is the background refractive 
index. 

The Maxwell-Bloch equations, together with the boundary 
conditions imposed by the laser cavity, provide the most 
general framework for studying laser instabilities. They are 
capable of handling both convective and absolute instabilities 
and can show transitions between them. However, their solu- 
tions require a numerical approach. Without host nonlinearities 
(7 = 0) and without GVD (/? 2 = 0), the steady-state solutions 
can be obtained, and their linear stability properties have been 
studied [9]. However, such an approach is quite cumbersome, 
and it is not easy to carry out the analysis after the inclusion 
of host nonlinearities and GVD. If one is interested only in 
absolute instabilities, an analytic approach can be developed, 
as discussed in the next section. 

III. Multimode Laser Equations 

Rather than solving (3)-(5) numerically, we make use of 
the fact that any cavity supports a set of longitudinal modes 
whose field distribution /m(*, 0 reproduces itself after each 
round trip inside the cavity. These modes can be obtained by 
solving (3) with Z? = 0 (no gain in the fiber cavity) and using 
the appropriate boundary irmditions at the cavity mirrors. For 
a high-Q laser cavity, one can distribute the mirror losses 
throughout the cavity and replace the fiber loss a in (3) with 
a? = a + QtAf, where (*m is the distributed mirror loss. The 
boundary condition then simply becomes A(L, t) — A(0, <)> 
where L is the cavity length. For a Fabry-Perot cavity with 
mirror reflectivities R x and ZZ 2 , a m is given by 



The approximation that the localized mirror loss can be 
replaced by a distributed loss only holds for a high-Q laser 
cavity [2]. 

When dealing with a unidirectional ring laser without host 
dispersion and nonlinearitics, the form of /m(*. t) becomes 
simply exp[— iw m (t — z/v g )\, where are the mode fre- 
quencies and the loss term has been ignored. In the presence of 
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GVD and SPM. / m (^, t) may depend on time in a complicated 
manner, but the temporal patterns remain invariant from round 
trip to round trip. However, the amplitude of such “temporal 
modes'’ can grow from round trip to round trip if an absolute 
instability occurs. It is thus important to consider time scales 
describing short-term and long-term dynamics, distinguished 
on the basis of the round-trip time T r = L/v g . The field- 
propagation equation (3) describes both the long-term and 
short-term effects. By making the coordinate transformation 
[10] 

r = <- 2 /t s , T = z/v g (7) 


mi ~Jc dr^ dT 

(13) 

/■Tr 


Cmnpq — I fmfnfpfq 

Jo 

(14) 


The values of these coefficients determine how many modes 
should be considered, i.e.. at what point the infinite sum in (8) 
and (9) can be truncated. 

The complex gain coefficient g in (11) is given by 

9mq — f g(T, T)f m {r)f^(r)dT (15) 


we distinguish between the propagation variable r describing 
evolution over a single round trip, and the long-term timescale 
T measured in units of T r . The steady-state solutions of (3)-(5) 
are found by putting the derivative with respect to T equal 
to zero and applying the boundary conditions at the cavity 
mirrors. 

In general, a laser can operate in several longitudinal modes 
simultaneously. Thus, the laser field can be written as A t {r) = 
Ylm If this steady state is perturbed the expansion 

coefficients a m become a function of T. We can study the 
long-term stability of the steady state by expanding the field 
.4 and the polarization B into a set of longitudinal modes 

<4(r.r) = £a m (r)/ m (r) (8) 

m 

B{T, r) =■ £ b m(T)f m (r). (9) 

m 


When the set of functions {/ m (r)} is complete, the modal 
expansion in (8) and (9) is exact. However, it is hard to prove 
completeness of { f Tn ( r) } under all operating conditions. As 
a practical matter, the infinite sum in (8) and (9) is always 
truncated. This is justified since only a finite number of 
longitudinal modes are excited in a laser at a given pump 
power. 

After expansions (8) and (9) are substituted in the 
Maxwell-Bloch equations (3) and (4) using the coordinate 
transformation (7). we multiply the equations for a m and b m 
by /*. and integrate them over one round-trip time T r . We 
then obtain 


1 v- dn m ^ 
2^ dT 


= \ Yi b ^ G ^q ~ \ “T ° 


m^mq 


->-E a m An, + h ^ &m C'mn pq ( 10 ) 


m, n p p 


m 

— (1 ^ ^ b rn G rng Z ^ ^ a mQrnq 


(11) 


where the coefficients G m9 . Z) mi7 . and C mnp? are defined as 



( 12 ) 


and should be interpreted as follows. When m — g, g mq 
represents the modal gain for the mode / m=? , while for m ^ q 
g mq represents gain modulation because of mode beating 
(sometimes referred to as population pulsations). From (5), 
we readily obtain the rate equation for the coefficients g m? : 

dT ~ ~~ 2 ~P~ — a n^m)C , m „p f . 

* m, n 

(16) 

Note that we have not made any assumptions about the 
orthogonal ity of the longitudinal modes. The nonorthogonality 
of longitudinal modes has attracted considerable attention 
in recent years [11]— [13]. In general, the more “open” a 
resonator is (the larger the mirror loss), the less orthogonal the 
longitudinal modes are. For fiber lasers, longitudinal modes arc 
expected to be nearly orthogonal since mirror losses are often 
kept to a minimum. 

If we assume that the functions {/ m } are orthonormal, the 
coefficients G mq reduce to 6 mq . The multimode equations 
(10), (11), and (16) then transform into the following set of 
Lorenz-Ha ken-type equations: 


1 da q i 1 i 

V t tr = 2 6 » - 2 QTa ’-2 /?2 E a - D ' 


mq 


+ *7 ^ a ™ a n a pC, 


mnpq 


m, n , p 


db 


T‘2 — — (1 — ib)bq — * y ^ dmffmg 


dT 

rp f 

ri -jjT —9oOpq “ Jpq 


~ “ < b m)C, 


mnpq • 


(17) 

(18) 

(19) 


Equations (17)-(19) form a complete set of equations de- 
scribing the dynamics of various longitudinal modes of the 
laser. They include the effects of host dispersion and nonlin- 
earity through the terms containing and j in (17). What 
constitutes SPM in a propagation-based description now splits 
into several different kinds of nonlinearities in the modal 
description. The triple sum in (10) describes the phenomenon 
of SPM when m = n = p = q, cross-phase modulation 
for m ~ q and n = p, and four- wave mixing for other 
combinations of m, n, and p [6]. 

Equations (17)-(19) constitute the main result of this paper. 
In principle, they can be used to study the phenomenon 
of passive mode locking by including a larger and larger 
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number of longitudinal modes as the mode -locked pulses 
become shorter. When the mode-locked pulse pattern becomes 
unstable, the laser may switch to a higher harmonic mode- 
locked operation [14] or even enter a chaotic regime [15]. 
Equations ( 17)— < 19) can be used to describe such behavior. 
When the steady-state pattern A $ (r) does not vary too rapidly, 
one can restrict the analysis to a small number of modes, 
although such a restriction will exclude passive mode -locking. 
We illustrate the usefulness of (1TH19) for this case by using 
a simple example. 

IV. A Simple Example 


where P -= X 2 is the (scaled) optical field intensity, Y is 
the (scaled) polarization amplitude. Aw, is the frequency shift 
with respect to the optical frequency wc, and <p t is the phase 
lag between electric field and polarization. We stress that all 
these quantities arc referring to the single longitudinal mode 
under consideration. 

After substituting (28)-(30) in the evolution equations 
(23}-{25) and denoting q = q f + iq*\ the field intensity 
P = X 2 is a positive real solution of the quartic equation 

4 

J2 c ” pn = 0 < 31 ) 

n=C 


The simplest case, analogous to the single-mode 
Lorenz-Haken model, is obtained from (17)-(19) by setting 
jji — n - p — q — 1. The resulting equations become 


1 da i 

Vg dT 

r dby 

T2 w 


* , «T , ,2 

-fc, - — ai + *7l a i I “l 

- (1 - »6)&i - I'aisn 


t 1 

Tl lr 


: <1 C - t/ll 


— Im(a*&i) 


( 20 ) 

( 21 ) 

( 22 ) 


where the SPM coefficient is defined as j — jC mi* The 
effect of GVD on a single longitudinal mode merely involves 
a frequency shift of magnitude w# = @ 2 V 9 Dn/ 2 ' which 
has been scaled out of the problem by changing the carrier 
frequency wd to wc + w#. This is equivalent to resetting the 
atomic detuning by 6 — <5+W£i. We normalize these equations 
using the standard Lorenz -Haken notation and rewrite them 
in the following way: 

a(r “ v) t- i<?x|x| 2 (23) 

(1 — ib)y -J- (r - z)x (24) 

bz -r Re(j**y) (25) 

where the normalized time t. — T jT> is measured in units of 
To. The parameters a. 6 , r, and q. and the variables x, y, and 
2 are related to quantities appearing in (3)— (5) as 


a — art-'jTb/S, 

b = Ts /Ti, 


»• - go/ oct, 

q - 7 P,TiV g , 

(26) 

x = y/b/P,a u 

V = (i&i/ar)\A/ Pt< 


2 =r - gu/oiT. 


(27) 


dr 

It 

dy 

dt 

dz 

dt 


These equations are identical to the standard Lorenz-Haken 
equations [3] except for the last term in (23) that is responsible 
for SPM and IDA occurring because of the host nonlinearities 

[3]. 

We now employ the standard linear stability analysis to 
investigate the stability properties of (23)-(25). First, we look 
for CW solutions of the form 


where the coefficients arc given by 


C 4 =q" 3 /<T 

(32) 

C 3 =[l+ 2 (l+ < r- l )k " 2 + (97^)kl 2 
c 2 =q"{tx+ Y) 2 /<x + 2 g"(! + a ~ l )( a + q") 

(33) 

- 2b6q'q"/a - brq" 2 + 6|g| 2 

(34) 

ci = [1 + (q" /<r)]{a + l ) 2 + b6 2 q" /a 


— 2bq"{(T + l)(r — 1 ) — 2b6q' 

(35) 

c c — bb 2 — &(r — 1 )(cr + l) 2 . 

(36) 

For each solution P , the accompanying frequency shift 
Aw,, inversion z t9 phase lag tp t9 and polarization Y arc 
determined by the following relations: 

. . 6 - q'P 

t.an cp, — Aw, + 5 = q»p 

(37) 

z, — r - f 1 + J (1 + tan 2 tp,) 

(38) 

v bz ‘ 

Y =~p coe <p,. 

(39) 


The presence of SPM affects the CW characteristics of 
the laser through the frequency shift Aw,. At resonance, i.c., 

5 = 0, the laser operates at the gain peak in the absence of 
SPM. In the presence of SPM, this CW state will be detuned 
from the gain peak, thus producing less output power. Thus, 
the effects of detuning 6 and the SPM parameter q f either 
counteract or re-enforce each other. In Fig. 1 we show how 
the laser power varies with pump at fixed SPM parameter 
q f = 1 for three different detunings 6. In general, the detuning 

6 increases the first laser threshold, but for 6 = 5 and r > 4 
the laser output power exceeds that obtained for 6 = 0. In Fig. 
1 , we choose a = 3 and 6 = 1 , typical values used in the 
context of the Lorenz-Haken model. 

The stability of the CW solutions is investigated by per- 
forming a linear stability analysis. After writing the perturbed 
steady states as 

x(t) = [X + «(<)] exp{-i[Aw,« + ip, + ^(0]} (40) 

lf(0 - \y + v(01 exp{-i[Aw,f + iMOl) ( 41 ) 

2(0 = + u>(<) (42) 


*,(0 =x exp[-i(Aw„« + <p,)] 
y f (t) - Y exp(— iAu/,0 
2,(0 = 2 , 


(28) 

(29) 

(30) 


linearizing the evolution equations for the perturbations 
u, v, u/, (where ^ is the phase difference Vu — and 

assuming a solution of the form exp(sf) for tx, t/, tt>, the 

characteristic eigenvalue equation for s is found to be a quartic 
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Fig. 1. The laser power P as a function of purnp level r in the presence of 
SPM (q* = 1) for three different values of atomic detuning 6 (indicated in 
the figure). Other parameters are a s 3 and b = 1. 

polynomial obtained from the determinant equation as shown 
in (43), shown at the bottom of the page. When any of the 
four complex roots of (43) has a positive real part, the steady 
stale under consideration is unstable. 

Because of the complexity of the determinant in (43), 
we do not attempt to find an analytical expression for the 
second threshold. Instead, we study two cases in detail. One 
is characterized by the combination <7 — 3 and 6 = 1, 
and is called the “bad-cavity” laser, while the other (“good- 
cavity” laser) has the combination cr — 6 = 1. The bad-cavity 
laser has a finite second threshold even in the absence of 
host noniinearities. while the good-cavity laser has no second 
threshold in the Lorenz-Haken limit obtained by setting 5 = 0 
[3j. For the resonant case ((5 = 0). the second threshold is 
simply given by r t h ,2 = cr(cr *+• b -h 3 )/( cr — 6 — 1) [3]. A 
necessary condition, the so-called bad-cavity condition, for the 
second threshold to exist is cr > 6+ 1. An analytical expression 
for the second threshold for arbitrary detuning values was 
given in [16]. 

We first consider a laser without IDA and set q n — 0. 
SPM is then the only host- induced nonlinearity. In Fig. 2, 
we show the pump parameter at the second threshold as a 
function of the SPM parameter q* for the bad-cavity laser 
operating on-resonance (b = 0). In absence of SPM ( q * — 0), 
the second threshold is located at r = 21. Quite surprisingly, 
the global effect of SPM is to stabilize the system: already at 
I?* I — 0.04 the second threshold ceases to exist. This result is 
surprising since SPM is essential for convective instabilities to 
occur in passive fibers and fiber amplifiers. From the steady - 



Fig. 2. The effect of SPM on the location of the second threshold at atomic 
resonance (6 = 0), for the “bad-cavity” laser with c = 3 and b = 1. 



Fig. 3. The effect of IDA on the location of the second threshold at atomic 
resonance (6 = 0). Solid line: bad-cavity laser with <7 = 1 and b = 1. Dashed 
line: good-cavity laser with <7 = 3 and b = 1. 

state analysis above we know that SPM causes a shift in 
operating frequency away from the gain peak, thus lowering 
the output power. From a stability point of view this shift can 
be interpreted as a shift away from the second threshold. The 
good-cavity laser (a = b = 1) shows no second threshold 
with or without SPM. 

Quite a different picture emerges when IDA is included. 
For both good- and bad-cavity lasers, EDA can reduce the 
second threshold dramatically, as shown in Fig. 3. In both 
cases, positive values of q ,f account for two-photon absorption 
while negative values of q H can be used to describe the 
effect of a fast saturable absorber. Clearly, fast saturable 
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Fig. 4. The location of the second threshold for the bad-cavity laser of Fig. 3, 
as a function of SPM for four different values of IDA (indicated in the figure). 

absorption (negative q H ) reduces the second laser threshold 
for both good- and bad-cavity lasers. In the presence of a fast 
saturable absorber, self-pulsing can therefore occur at pump 
levels only a few times above the first threshold. We stress 
that this self-pulsing is not related to passive mode locking 
because of our assumption of a single longitudinal mode. In 
fact, the frequency 7 of self-pulsing is related to the relaxation 
oscillations and is a fraction of the longitudinal-mode spacing. 

Fig. 4 shows the combined effect of SPM and IDA for the 
bad-cavity laser. For four different values of q*\ the second 
laser threshold is shown as a function of the SPM parameter 
q* . The effects of SPM and negative IDA have opposite effects 
on the second threshold: SPM stabilizes the laser by increasing 
the second threshold, while negative EDA destabilizes the laser 
and reduces the pumping level at which the second threshold 
occurs. For the good-cavity laser, a qualitatively similar figure 
is obtained, although in that case no second threshold exists 
when q " = 0 . 

Fiber ring lasers usually oscillate in many longitudinal 
modes simultaneously because their mode spacing (~10 MHz) 
is a fraction of the gain bandwidth (>1 THz). By using a 
grating, they can be forced to operate in a single longitudinal 
mode. When we use typical values for a neodymium-doped 
fiber laser, aj — 46 km -1 , Tn — 1.75 ps, T\ — 0.1 ms, 
7 — 2.3 W" l km“ l , and P t — 11.4 /iW, the dimensionless 
parameters are a = 8.84 x 10 -6 , b — 1.77 x 10“ 8 , and q f — 
5.23 x 10” 4 . These values make the single-longitudinal-mode 
fiber ring laser a good-cavity laser with no second threshold by 
itself. Fig. 5 shows the effect of fast saturable absorption on the 
second threshold of such a laser. Even relatively small values 
of q n cause a second threshold to exist at quite low pump 
values. For example, when q" — —3.65 x 10~ 8 , the second 
laser threshold exists at a pump level only twice the first laser 
threshold. This value of q H corresponds to Im( 7 ) = 1.6 X 10 ~ 4 
W^km " 1 if we use C mi — 1. This is a weak saturable 
absorption indeed, since it only reduces the losses in the cavity 
by 0.001 % for an intracavity power of 3 W. At the second laser 
threshold, the relaxation-oscillation frequency 7 is 50 kHz, and 
self-pulsing will occur with a repetition rate close to that value. 



Fig. 5. The location of the second threshold for a fiber ring laser which is 
forced to oscillate in only one longitudinal inode, as a function of (negative) 
I DA. 

V. Conclusions 

In this paper, we have analyzed the effects of GVD and 
host nonlinearities (SPM and IDA) on the absolute stability 
of lasers. Starting from the Maxwell-Bloch equations, we 
derive a set of multimode laser equations, which govern the 
absolute dynamics of the laser on a timescale longer than 
the cavity round-trip time. The presence of a resonator gives 
rise to longitudinal modes that can differ from their conven- 
tional form in the presence of GVD and host-nonlinearitics. 
Our analysis shows the link between the propagation-based 
(convective) modulation instability 7 and the purely temporal 
(absolute) instabilities of these spatio-temporal patterns. The 
latter can occur only in a resonator as absolute instabilities 
require optical feedback. 

To illustrate our formalism, we consider the case of a fiber 
ring laser forced to oscillate in a single longitudinal mode. The 
nonlinear dynamics of the slowly varying amplitude of that 
mode is governed by the Lorenz-Haken equations except that 
a nonlinear term is added to the field equation. This nonlinear 
term contains a complex parameter q whose real and imaginary 
parts account for nonlinear effects such as SPM and IDA, 
respectively. SPM is found to stabilize the laser, while IDA, 
in general, destabilizes the4aser. Thus, we show that SPM can 
have a double role as far as stability issues are concerned. 
On the one hand, it can cause convective instabilities in 
combination with GVD, which can lead to exponential growth 
of a localized perturbation. However, once the perturbation 
fills the entire cavity, its subsequent growth is dictated by 
the cavity modes. SPM is found to suppress such subsequent 
growth. EDA induced by fast saturable absorbers, however, 
quickly gives rise to a second laser threshold at relatively 
low pump levels. In our single-mode laser, the laser begins to 
self -pulse at a repetition rate close to the relaxation-oscillation 
frequency. 
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